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Abstract. After showing that the anelastic approximation is better 
than the subseismic one to filter out acoustic waves when studying low- 
frequency stellar oscillations, we compute gravito-inertial modes of a typ- 
ical 7-Doradus star using this approximation. We show that eigenmodes 
can be regular or singular, according to the possible focusing towards 
attractors of the underlying characteristics. Consequences on the oscilla- 
tions spectrum are then discussed. 



1. Introduction 

7-Doradus stars are rapidly rotating variables for which the second-order per- 
turbative theory of Ledoux (1951) fails to reproduce observations. Indeed, os- 
cillations and rotation periods are both around one day (Handler & Krisciunas 
1997) and eigenfrequencies should be non-perturbatively computed. 

In this paper, gravito-inertial modes of a typical 7-Dor star are computed 
using the anelastic approximation. We first compare the anelastic and subseis- 
mic approximations and show that the anelastic one is the best to filter out 
acoustic waves from the infinite system of differential equations (§2). Properties 
of gravito-inertial oscillations are then recalled (§3) and the part played by the 
characteristics trajectories is enlightened, before conclusions in §4. 

2. Comparison of the anelastic and subseismic approximations 

These two approximations have been compared in Dintrans & Rieutord (2001). 
Both approximations neglect perturbations in the gravitational potential and 
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assume that the Eulerian pressure fluctuations do not contribute to Lagrangian 
ones. On the contrary, they differ on the form of the equation of mass conser- 
vation as 

Anelastic: • (pov) = 0, Subseismic: ■ v = {g/(?)vr, 

where p^, g and respectively denote the equihbrium density, gravity and 
square of the velocity of sound whereas v is the velocity. 

Comparisons have been made using g-modes of the homogeneous and n = 3 
polytropcs. Analytic expressions for the anelastic and subseismic eigenfrequcn- 
cies have been found for the homogeneous model whereas all the eigenfrequencies 
have been computed numerically for the n = 3 polytrope. 

Results are summarized in Table 1. They show that the anelastic values 
are between 5 and 20 times more accurate than the corresponding subseismic 
ones demonstrating that the anelastic approximation is best to filter out acoustic 
waves when studying low-frequency stellar oscillations. 



Table 1. Complete dimensionlcss eigenfrequencies lS^ for some g- 
modes of the = 3 and homogeneous (in parenthesis) polytropes, with 
their associated anelastic w^^gj and subseismic Wg^^ counterparts. 





X 10^ 


^ancl X 


^sub X 


g5 

gio 

g20 
g30 


5.76059 (-36.1367) 
1.98396 (-12.2196) 
0.59984 (-3.65629) 
0.28631 (-1.73527) 


6.03954 (-35.2941) 
2.02646 (-12.1212) 
0.60453 (-3.64741) 
0.28744 (-1.73327) 


6.27794 (-32.6087) 
2.08345 (-11.5607) 
0.61587 (-3.55239) 
0.29154 (-1.70180) 



3. Application to the long-period oscillations of a 7-Dor type-star 

The anelastic equations need to be solved in a co-rotating frame; they read 
(Dintrans & Rieutord 2000, hereafter DR2000) 



ia'^ jv ■erj+'^x (2lt x i; ) = and ^ • {pqv) = 0, (1) 

where we assumed a time-dependence of the form eyq){iat) and N denotes the 
Brunt- Vaisala frequency. After projecting on vectorial spherical harmonics (see 
Dintrans, Rieutord & Valdettaro 1999 for numerical details), one obtains a gen- 
eralized eigenvalue problem which is solved using either the QZ algorithm (all 
of eigenvalues a arc computed) or an iterative Arnoldi-Chcbyshcv solver (only 
desired eigenvalues (7„ are computed, with their associated eigenvectors Vn)- 

The perturbative theory has been tested by following a high-order g-mode 
(pcriod'-^0.3d) of a typical 1.5 Mq 7-Dor star as rotation was increased. As 
expected for this kind of long-period mode, perturbative theories rapidly reach 
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Figure 1. Integration of the characteristics equation (^) for a star 
rotation period of one day and two mode frequencies a. If a is well 
above 2Q (a), characteristics fill the whole hyperbolic domain whereas 
they can be focused along attractors when a ~ 2fl (b). Three hundred 
reflections have been drawn in both cases but we removed the first two 
hundred ones in plot (b) to emphasize the final attractor. 

their limits since first- and second-order relative errors are comparable to the 
observational detection limit for rotation periods of about 8 and 3 days, respec- 
tively. 

Then, by solving the non-perturbative problem (jl]) in the rapid rotation 
regime (rotation down to one day), we found the more and more important part 
played by the characteristics trajectories as the rotation increases. Characteris- 
tics, on which energy propagates, obey the following differential equation (s and 
z being cylindrical coordinates with = + z^) 

dz _ N^sz ± r^a'^N'^s'^ + (4J^^ - a'^){a'^r^ - N'^z^) 

Given a mode frequency a and rotation rate Jl, this equation can be integrated 
in the radiative zone of the star and we found that the resulting web imposes 
the mode structure depending on the ratio a/ 20.: 

• if (7 ^ 2ri, all characteristics trajectories are ergodic (Fig. la) and the 
associated eigenmodes are regular, in the sense that the velocity field is 
smooth and square- integrable (see Fig. 8a in DR2000). 

• if o" ~ 2Vt (the regime of interest), characteristics trajectories almost focus 
along attractors (Fig. lb) and the associated eigenmodes are singular, that 
is, velocity tends to infinity on the attractor (see Fig. 8b in DR2000). Thus, 
without an efficient diffusion process, these modes cannot exist. 

• if o" < 2r2, no global eigenmodes have been found, most likely because of 
numerical difficulties. 
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4. Conclusion 

First, we showed that the anelastic approximation is better than the subseismic 
one to filter out acoustic waves from oscillations equations. Second, wc applied 
this result to investigate the low-frequency oscillations of a rapidly rotating 7- 
Dor type-star. 

We tested first- and second-order perturbative theories by following a long- 
period g-mode while increasing rotation. We showed that a non-perturbative 
approach is necessary for rotation periods below 3 days. Then, a computation 
of eigenmodes allowed us to reach the one-day rotation regime where we found 
that modes part into regular and singular ones as the associated characteristic 
trajectories are ergodic or focused along an attractor. It has important conse- 
quences on the oscillations spectrum since frequency subintervals corresponding 
to attractors are lacking eigenvalues when dissipative effects are neglected. 
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Discussion 

R. Townsend : What advantage does the anelastic approximation give, if you 
must still use an expansion in spherical harmonics ? 

B. Dintrans : Wc always need to expand the angular part of perturbations onto 
spherical harmonics. But adding rotation through Coriolis' force leads to an 
infinite system of coupled differential equations. The usefulness of the anelastic 
approximation is precisely to decrease the system size by filtering out acoustic 
waves, which makes numerics more efficient. 

R. Townsend : What is the physical interpretation of the singular modes which 
you find ? 

B. Dintrans : These modes are due to the focusing of the characteristics, as- 
sociated with the inviscid problem, towards periodic orbits (attractors); this 
process is therefore coming from the ill-posedness of the associated mathemat- 
ical problem. Physically, you may think to a wave packet which, following the 
characteristics, will bounce on the boundaries. After each cycle, its wavelength 
is reduced by a factor (which depends on frequency) larger than unity. After an 
infinite time, the wavelength is zero and the singularity is born. 



